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The worldsheet formulation is introduced for lattice gauge theories with dy- 
namical fermions. The partition function of lattice compact QED with staggered 
fermions is expressed as a sum over surfaces with border on self-avoiding fermionic 
paths. The surfaces correspond to the world sheets of loop-like pure electric flux 
excitations and meson-like configurations (open electric fiux tubes carrying matter 
fields at their ends). The proposed formulation does not have the problem of the 
additional doubling of the fermion species due to the discretization of time. The 
gauge non-redundancy and the geometric transparency are two appealing features 
of this description. From the computational point of view, the partition function 
involves fewer degrees of freedom than the Kogut-Susskind formulation and offers 
an alternative and more economic framework to perform numerical computations 
with dynamical fermions. 
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I. INTRODUCTION 



The problem of handling dynamical fermions is still a major challenge that faces lattice gauge 
theory at present. Monte Carlo techniques have provided many important results clarifying several 
points about the dynamic, and recent computations are achieving 10 % or better accuracy in the 
spectrum both for heavy quark and light quark system [0. However, statistical algorithms are 
very expensive in computer time, and the increased computing power of the coming generation 
of machines without further theoretical insight will probably be insufficient in order to definitely 
improve results. Basically, the difficulty posed by the fermions stems from the fact that they are 
represented not by ordinary numbers but by anti-commuting Grassmann numbers which cannot 
be directly simulated numerically. Since the fermion field appears quadratically in the action, the 
usual procedure is to integrate it out producing the Mattews-Salam determinant. So, the problem 
of including dynamical fermions is reduced to one of evaluating the determinant of a large matrix. 
This is a costly task. 

An alternative to tackle this problem is to resort to the analytical methods. These could be 
divided into two categories: strong-coupling expansion and Hamiltonian variational-like methods. 
The principal limitation of the former is the difficulty to reach the weak-coupling region. An 
exponent of the second group is the loop approach ||^. The basis of the loop method can be 
traced to the idea of describing gauge theories explicitly in terms of Wilson loops or holonomies jj]- 
since Yang |^ noticed their important role for a complete description of gauge theories. The loop 
Hamiltonian was given in terms of two fundamental operators: the Wilson Loop operator (the trace 
of the holonomy) and the electric field operator (its temporal loop derivative). The loops replace 
the information furnished by the vector potential (the connection) . A description in terms of loops 
or strings, besides the general advantage of only involving the gauge invariant physical excitations, 
is appealing because all the gauge invariant operators have a simple geometrical meaning when 
realized in the loop space. Last but not least, the interest on loops relies on the fact that it was 
realized that this formalism goes beyond a simple gauge invariant description and in fact it provides 
a natural geometrical framework to treat gauge theories and quantum gravity. The introduction 
by Ashtekar Q of a new set of variables that cast general relativity in the same language as gauge 
theories allowed to apply loop techniques as a natural non-perturbative description of Einstein's 
theory. 

In 1991 the loop representation was extended in such a way to include dynamical staggered [^j 
fermions : the so-called P-representation Roughly the idea is to add to the closed pure gauge 
excitations, open ones corresponding to "electromesons" . Afterwards the P-representation was used 
to perform analytical Hamiltonian calculations, by means of a cluster approximation, providing 
qualitatively good results for the (2-1-1) |10| and the (3 -I- 1) [|ll| cases when compared with the 
standard Lagrangian numerical simulations in terms of the fields. The Hamiltonian method has 
the serious drawback of the explosive proliferation of clusters with the order of the approximation. 

Thus, our goal was to explore another approach: to build a classical action in terms of strings 
and knit together the transparency and non-redundance of the string P-formulation and the power 
of the Lagrangian simulations. The first step of this program was the introduction of new lattice 
action for pure QED in terms of closed strings of electric flux (loops) |12|. In the pure case the 
action is written as a sum of integer variables attached to the closed worldsheets of the loop 
excitations. The second step was to include matter fields into the string description, with this aim 
we considered the simplest gauge theory: the compact scalar electrodynamics (SQED) In the 
case of SQED the action is expressed in terms of open and closed surfaces which correspond to 
world sheets of loop-like pure electric fiux excitations and open electric flux tubes carrying matter 
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fields at their ends. The previous two worldsheets actions were simulated using the Metropolis 
algorithm being the results quite encouraging , Q . 

Here we show how to introduce dynamical fermions in a worldsheet or Lagrangian description. 
This paper is organized as follows. Secction II is a 'bird-eye- view' review of the P-representation on 
a hypercubic lattice. We show the realization of the lattice QED Hamiltonian in the Hilbert space 
of paths {P}- In section III we present the Lagrangian counterpart of the previous P-representation 
and we write the partition function in terms of worldsheets of the string-like excitations. In section 
IV, by using the transfer matrix procedure, we check that we get the Hamiltonian of section II from 
the path integral of section III. Notice that this enabled us to get the Kogut-Susskind formulation 
via the transfer matrix, this is an interesting problem which was not properly solved. Finally, 
section V is devoted to conclusions and some remarks. 



II. THE P-REPRESENTATION ON THE LATTICE 

The P-representation offers a gauge invariant description of physical states in terms of kets 
I P >, where P labels a set of connected paths PJ' with ends x and Q In order to make the 
connection on the lattice between the P-representation and the ordinary representation, in terms 
of the fermion fields ip and the gauge fields C/^(a;) = exp[ieaAp(a;)], we need a gauge invariant 
object constructed from them. The most natural candidate in the continuum is 

<^>{py)^^p\x)u{py)4'{y), (i) 

where U{py) — exp[jea JpA^dx^]. 

The immediate problem we face is that $ is not purely an object belonging to the "configu- 
ration" basis because it includes the canonical conjugate momentum of "0, ij}'^ . The lattice offers 
a solution to this problem consisting in the decomposition of the fermionic degrees of freedom. 
Let us consider the Hilbert space of kets | tp'ljtpd, A/^i >, where u corresponds to the up part of 
the Dirac spinor and d to the down part. Those kets are well defined in terms of "configuration" 
variables (the canonical conjugate momenta of tpd and i/'J are i/'d 4'u respectively. Then, the 
internal product of one of such kets with one of the path dependent representation (characterized 
by a lattice path PJ' with ends x and y) is given by 

^py)^<py-:n^l^d,A^> 
= €A^)uiP^)^M (2) 

where i and j denote a component of the spinor u and d respectively. Thus, it seems that the choice 
of staggered fermions is the natural one in order to build the lattice P-representation. Therefore, 
the lattice paths PJ* start in sites a: of a given parity and end in sites y with opposite parity. The 



^For a more detailed exposition of the P-representation and the realization of the different operators see 
the ref. §. 
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one spinor component at each site can be described in terms of the xi^) single Grassmann fields 
The path creation operator <i>Q in the space of kets {| P >} of a path with ends x and y is 
defined as 

^Q=xH^)UiQ}^my)- (3) 

Its adjoint operator <i>g acts in two possible ways annhilating the path or joining two 
existing paths in | P > one ending at x and the other strarting at y. 

Let us show the realization of the QED Hamiltonian in the Hilbert space of kets \P >■ This 
Hamiltonian is given by 

H = {g^l2)W 

W = We + XWi + \^Wm 
X 1/g^ 

e 

Wj = -^($, + $1) 

i 

($1 = Tjn{x)xH^)Un{x)x{x + n)) (4) 
??-e.(a; + aei) ^ Veiix), 
p 

where x labels sites, £ = {x, n) the spatial links pointing along the spatial unit vector n, 
p = (x,n,n') the spatial plaquettes; is the electric field operator, which is diagonal in the 
P-representation 

Ei\P >= Ni{P)\P >, (5) 

where the eigenvalues Ni{P) are the number of times that the link I appears in the set of paths 
-P; Up = Wi^pUi- The $f are "displacement" operators corresponding to the quantity defined in 
for the case of a one-link path i.e. P = £. The realization of the different Hamiltonian terms 
in this representation is as follows 

First, by (||) the action of the electric Hamiltonian is given by 

We I P>^J2 ^«^(^) I ^ > • 
The interaction term Wj can be written as 
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where the subscripts e and o denote the parity of the lattice sites. This term is reaUzed in {| P >} 
as 

-Wj\P>^Y^ e(P,4J I P • 4. > + 51 <^'>4J I P • 4„ > (7) 

where 4 is the Hnk starting in x and ending in x + na. For Hnks of even origin, e{P,£) is zero 
whenever an end of £ coincides with an end of P, it is —1 when $£ "deletes" the link i ^ P , 
dividing one connected component into two, and it is +1 in any other case. For links of odd origin, 
e(P, £) is zero unless both ends of £ coincide with two ends in P. In that case, it is —1 when £ 
joins two disconnected pieces and it is +1 when it closes a connected piece or when it annhilates a 
link. The different actions of operators over path-states \P{t) > are schematically summarized 
in FIG.l. 




(a) (b) (c) (d) 

FIG. 1. A summary of the different actions of operators "1>£ applied over path-states \P >. The link £ is 
represented by a dashed bond in the pictures on the top and the resulting paths \P' > are plotted below. 

Finally Wm is the sum of the operators U{p) and U^{p) which add plaquettes, and can be 
written asQ 

p 

where 



^In such a way that a generic path P is generated out from the 0-path state (strong coupling vacuum) 
by the application of the operator string xU^) n^GP" VeUix{y). 
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and then 

Wm\P>=Y.^\P-P> + \P-P>) (8) 
p 

where p and p respectively denote the clockwise and the counter-clockwise plaquette contour. 



III. THE WORLDSHEET OR P-ACTION 



In order to cast the preceding path description in the Lagrangian formalism let us begin by 
considering the path integral for lattice QED with staggered fermions. We will show that this 
leads to a surface action which is not possible to connect directly with the Hamiltonian via the 
transfer matrix in the Hilbert space {\P >}. However, this action will serve as a guide in order to 
build the genuine P-action Sp. 

For simplicity we choose the Villain form of the action, which is given by 



Z= jw^dx] /'[d^]^exp{-^^[0p + 2™p]2 + l^a(^-i)ry,(x];C/«.+aA + ^-C.)} 

= jw^dx] l[de]Y^e^^{-^\\ve + 2^nf + ]-Y,a^''-^^m{xlUiXr+a[. + H.c.)}, (9) 

where we used in the second line the notations of the calculus of differential forms on the lattice 
of [Q. In the above expression: D is the lattice dimension, /3 = ^ is a real compact 1-form 
defined in each link of the lattice, Ui = e*^'^ and % and x^ ^-re Grassmannian variables defined on 
the sites of the lattice, V is the co-boundary operator -which maps A:-forms into (fc -I- l)-forms -, 
n are integer 2-forms defined at the lattice plaquettes and =< g,g >= J2ck 5^(cfc); where g 
is any /c-form and are the /c-cells (cq sites, ci links, . . . ). The measure in (j^) is 

[dx^dx] = n "^^rdXr, m = n 5 (10) 

r e 

Let us forget for the moment the global factor a^^^^K The equation (^) can be written as 

Z= j[d9]cM-Sgauae{O)}ZFi0), (11) 

where 

Zf^ f [dx^dx] exp{i ^?7^(x];;7,Xr+aA + H.C.)}, (12) 

where Sgauge stands for the pure gauge part and Zp denotes the fcrmionic partition function. Now 
if we expand the exponential in Zp we get 
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Zf^ [dx^dx]'[[exp{-T]exlU(:Xr+aii}Y[exp{-r]exl+af,ulxr} 

= f [dx^dx] Hii + \mxluixr+ai.) n(i + \mxUa-,ulxr). (13) 

i.e. in the above product we have to consider each link and its opposite. 
Let us recall the rules of Grassmann variables calculus: 

0, (14) 

0, (15) 

1. (16) 

Therefore, when we expand the products in (|l^) and the Grassmann variables are integrated 
out, the only non- vanishing contributions arise from these terms with x\ and Xr appearing one and 
only one time for every site r. In other words, the integration of the Grassmann variables produce 
products of Uis along closed paths i.e. Wilson loops. We denote by a generic configuration of 
multicomponent paths. T is specified by a set of oriented links verifying the rule that they enter 
and leave one and only one time each lattice site r. This self-avoiding character is the geometric 
expression of the Pauli exclusion principle. 

We distinguish two parts in !F: the set of connected closed paths where a link is never run 
in more than once, JF'^ and the set of isolated links traversed in both opposite directions or "null" 
links T. Then T ~ T"^ \J T . The number of connected closed components of each kind is called 
Njrc and N jp. In FIG. 2 we show a two-dimensional sketch on a 2 x 2 x 3 lattice of a possible 
configuration T consisting in one fermionic loop and two "null" links. 



{Xv,Xs] = {Xr,xl} = {Xr,xl} 
dXr = \ dxl 



dXrXr = / dxlx 



A- 

1 1 






r - 







I I 

FIG. 2. A possible configuration of self-avoiding paths in a 2 x 2 x 3 lattice. Fermionic loops are 
represented by filled bold lines and "null" links by dashed bold lines. 



The Grassmann integration over each element of T'^ or J- gives a —1. Thus, after integrating 
the fermion fields, the fermionic path integral becomes 



z^=^(-i)^-(-i)^^n'?^f^^' 



(17) 
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where V is the total number of lattice sites and we used that the terms rjeUe for £ in the part J- 
of T cancel out. It is easy to check that the number of "null" links in a given T is connected 
with the number Nyrc of fermionic loops as follows^ 



= (18) 



where Ljra is the number of links in T'^ (the length). Thus, up to a global sign, we have 

(-1)^^ = (-1)"^. (19) 
One can define the ?7jrc for a fermionic loop T'^: 

rir^^Wm^ n r;p = (-l)^-% (20) 
leJ^" peS:Fc 

where Aj^c is the number of plaquettes which make up any surface Sj^i enclosed by JF'^ (different 
choices of SyrcSj^c differ by a even number of plaquettes so the sign is well defined and one can 
choose the Aj^c as the minimal area bounded by JT'^ ) . Neglecting a global sign we get 

c 

ZF = J2<Jr'[[U,, (21) 

where 

= (-l)^^'^-^+^^^ (22) 

Let us analyze a little closer the sign ajr. For D — 2 we show in the APPENDIX I that all the 
non- vanishing contributions have cr^r = +1. In more than two space-time dimensions there are 
non-null contributions J- to Zp with both signs. The reason is that connected fermionic loops 
enclosing odd numbers of vertices do not necessarily imply any more a null contribution. Different 
examples of a simple fermionic loop contributing with For instance, in case (a) the enclosed area 
Ajrc is 4, the lenght Ljc-c is 8. Therefore, the fermionic loop (a) has TVjrc — -f Ajrc = 1 — 4-1-4 = 1 
, and then ajr = — 1. 



^This is true for a lattice with an even number of sites (a lattice with an odd number of sites can not be 
"filled" with T paths). 
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(a) (b) (c) 

FIG. 3. Two examples of fermionic loops contributing with a -1. 
Coining back to the total path integral we can write it as: 

We can express the fermionic paths T'^ in terms of integer 1-forms -attached to the hnks- / 
with three possible values: and ± 1 with the constraint that they are non self-crossing and closed 
9/ = where d is the boundary operator adjoint of V which maps fc-forms into (fc — l)-forms. 
Both operators verify the integration by parts rule 

< 5, V/i >, (24) 



where g and h are respectively fc and fc — 1 arbitrary forms. The factor Jl/ej^<: ( P3[ ) is nothing 
but the product of Wilson loops along the f- loops i.e. exp{i 'YliteT" ^^/^}- ^ can be expressed 
by means of functions attached to links / with value or 1 and defined over /^^(O), the "nucleus" 
of /. In terms of the / and / we get: 

^= /Me]EEE°^P{-fll^^ + 2™f + *<^'/>}- (25) 

" / / 

If we use the Poisson summation formula g(n) — dBg{B)e^'^^^^ -where n and n' are 

integer 2-forms and i? is a real 2-form- we get 

z= /['^^]EEE^(/)x 
" / / 

/ [dB]exp{-^||V6l + 27rBf + i < 6*,/ > +27ni < B,n >}. (26) 
Performing the displacement B ^ B — V9/2tt and integrating in B 

Z = I m '^(/) exp{-^||nf + ^<eJ-^n>} 

^ f f 



^(/) exp{-^||nf }'5(/ - dn) 

^ f f 



= ^X^c7(5n)exp{-^|lnf}, (27) 

n ~f 

where we have used equation ( p4[ ) in order to transform < V0, n > into < 0, dn > whose integration 
produced the Dirac's delta 6{f — dn). Equation ( p7| ) is a geometrical expression of the path integral 
of lattice QED with staggered fermions in terms of surfaces with self-avoiding boundaries. 

Now, let us return to our goal, namely to set up the worldsheet Lagrangian formulation 
corresponding to the P-representation of sect. 2. The path integral ( p7| ) includes worldsheets of 
paths with ends of any parity. Therefore, obviously, it does not produce the Hamiltonian (|^) by 
the transfer matrix method. 

As a matter of fact, the transfer matrix for staggered fermions is a general and interesting 
problem on which not much has been done. Unlike the Wilson fermions which come out rather 



nicely from the transfer matrix ||l5[- |17|, the staggered fermions present some troubles. For the 
spatial part of the action the derivation goes in parallel with that for the Wilson theory. The 
interesting question is how does the temporal part work. This problem was discussed in the work 
of Saratchandra et all JTst , in this paper was analyzed the transfer matrix formalism for relating the 
Kogut-Susskind Hamiltonian and the Euclidean action for staggered fermion fields. The authors 
pointed out that the ordinary case, in which the one-component fermions xi^) ^^id X^{^) both live 
on all sites, exhibits a doubling of the fermion species with respect to the hamiltonian formalism 
(four flavours instead of two flavours). Additionally, the corresponding transfer matrix is hermitean 
but not positive definite. They explored the alternative of considering x and on alternate sites, 
although this procedure avoids the excess of flavours and gives a positive transfer matrix it is not 
free from complications. 

Although the partition function (27) is not directly connected with the Hamiltonian (^ via 
the transfer matrix, it is closely related with it (as we will show) and provides us with a guide to 
guess the genuine partition function for the P-representation Zp. With this aim we consider the 
restriction on the set of the surfaces, to the subset of the surfaces such that they are worldsheets 
of the P paths. This is equivalent to require that when intersected with a time t — constant plane 
they give the paths of {\P >} i.e. paths with ends of opposite parity and oriented from even sites 
to odd sites. We get a link of Pt for every plaquette of the surface which connects the slice t with 
the t + ag. As an illustration, in FIG. 4 we show the paths of the P-representation we get from the 
minimal surface enclosed by the configuration of self-avoiding T'^ paths depicted in FIG. 2. 



I 

J -> 

' I ' I 

I I I I 

I 



I 
I 
I 



>r 



FIG. 4. The P-paths (thick lines) for the different times corresponding to the surfaces enclosed by the 
configuration of self-avoiding paths of FIG. 2. 

The imposed constraint forbids fermionic loops such that the one of FIG.3-(c) which produces at 
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the first temporal slice a path connecting vertices with the same parity. The first configuration 
with negative Boltzman factor which appears in the partition function, i.e. that of lower action, is 
the spatial square of side 2 depicted in FIG.3-(a) with area Ajrc = 4. 
Thus, we propose the following worldsheet partition function: 

Zp = J2MS)eM-^Y.<}^ (28) 

s pes 

where S runs over worldsheets of P-paths and ah{S) is a sign defined: 

'^h{S) - i^iPt-idStHPt) (29) 

t 

where St are the surfaces given by the spatial plaquettes of S at section t, At is the number of 
plaquettes in St and the function gives a sign defined by 

HP) = n(-i)^ n(-i)^ (30) 

a p 

[P — {PcnQp} being P^ and the open and closed component paths respectively. See AP- 
PENDIX II for notation and discussion.) Additionally, note that we have eliminated the sum 
over the "null" links. The rationale for this is that the "null" links do not play any role in the 
P-representation. In the next section we are going to show explicitly that Zp gives rise to the 
Hamiltonian (^) via the transfer matrix procedure. 



IV. HAMILTONIAN OBTAINED VIA THE TRANSFER MATRIX METHOD 



By means of the transfer matrix method let us show that we re-obtain the Hamiltonian (|j) 
from the path integral Zp. As we wish to consider the continuous time limit of the previous lattice 
Euclidean space-time theory, we introduce a different lattice spacing oq for the time direction. The 
couplings on timelike and spacelike plaquettes are no longer equal in the action i.e. we have two 
coupling constants: /3o and (3s- The temporal coupling constant (3q decreases with ao whilst the 

spatial coupling constant (3s increases with oq. So far we have neglected a factor — in the 
path integral. Taking into account the fact that the lattice has a different temporal separation 
ao ^ a we get a relative factor of (^)'''^'' for each spatial link of T'^. 

s ^ pes 

where \dn\sp denotes the number of spatial links in dn. To factorize Zp to fixed time contributions 
we consider the spatial plaquettes of S that for each t define the spatial surfaces St and the temporal 
plaquettes that define the spatial paths Pt and we write 
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and 

'Jh{S) = l[<7'{Pt-,,Pt,St) (33) 
t 

a'{Pt^i,PuSt) = {-l)^'j^{Pt-idStHPt) (34) 

To write the operator T which connects the ket \Pt-i > with the ket \Pt > we begin by 
decomposing the sum over world sheets J^s (HI) ^ sums: one over the temporal parts, the 
Pt, and one over spatial parts St i.e. 

E n< pt\T\Pt-i > 
{Pt} t 

< P'\f\ P>^Y1 «r'"'''^V*exp{--L||pf - } (35) 

Taking into acount the relation which connects T with H when ao is small: 

f oc e-'^«^+°('^o\ (36) 

we find H considering this limit 

< P'\f\ P >« Sp^p, -ao< P'\H\ P > 

The dominant contributions in this limit are the cases where P' is equal to P or differs of it by 
one link £ or by one plaquette p. The power of ao in ( ^5| ) forces St ^ in the first two cases and 
5f = p in the third. 

< P|f|P>«l«exp{--l-||Pf} (37) 

^po 

< P ■ p\f\ P>^-ao< P- p\H\ P >« - exp{--^} (38) 

< P ■ e\f\P >K -ao < P ■ e\H\P >^ aoa* (39) 
To obtain a proper continuum time limit we should take 

/3o = ^ (40) 
" 2ln{2g^a/ao) ' ^^^^ 



12 



where a continues to denote the spacehke spacing. The values of cr* are — 1 when adding a plaquette 
( |38| ) and depend on the way we do it when adding a hnk (^9|). An anahsis of the possible cases 
shows that cr' corresponds to the sign e in (^. 
Then, 



SO we recover the Hamiltonian (jj) . confirming that (^8|) is the expression of the partition function 
of compact electrodynamics in the P-representation. 

In FIG. 5 we show an scheme which summarizes the different actions/hamiltonians for QED 
with staggered fermions, and their connections. The surface action >5'gyj.-f of (^7|), which was 
obtained from the Kogut-Susskind action Sks by integrating over the gauge and fermion fields, 
does not produce the hamiltonian Hp which in fact is connected with the Sp of (p8|). This is 
because Sp is free from the known problem of Sks of having an unwanted additional doubling 
compared with the hamiltonian formulation. On the other hand, the equivalence between the 
Kogut-Susskind Hamiltonian Hks and the P-Hamiltonian Hp was proved in |9|. 




FIG. 5. A diagram summarizing the different lattice actions and Hamiltonians for compact QED and 
their connections. The horizontal lines denote equivalence. The vertical lines emanating from the classical 
actions to the hamiltonians represent the connection via the transfer matrix procedure. 

In other words, the action Sp should be regarded as a different lattice action, which produce 
a Hamiltonian equivalent to the Kogut-Susskind's one. 



V. CONCLUSIONS AND FINAL REMARKS 



We propose a purely geometric action in terms of the world-sheets of P-path configurations. In 
this formulation, the connection between the Hamiltonian and Lagrangian of QED with staggered 
fermions is straigthforward via the transfer matrix method. The partition function can be written 
as a sum over surfaces with border on fermionic self-avoiding loops . Hence, the fermionic 
problem has been reduced to the task of computing quadratic areas enclosed by polymer-like 
configurations. The polymer representation of lattice fermions |ig| ] is often used in a different way 
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to compute the fermionic determinant. It is important to note that our formahsm is free from the 
problem of the additional doubling of the fermion species due to the discretization of time. 

With regard to the economy and possible advantages from the numerical computation point 
of view offered by the P-description, we want to emphasize two facts: I) Concerning the gauge 
degrees of freedom, it only involves sums of gauge invariant variables i.e. no gauge redundancy. 
II) It includes a subset of the configurations which are taken into account in the path integral ( P7|) 
equivalent to the one of Kogut-Susskind. Our formulation involves a sum of configurations with 
Boltzman factors of both signs. However, the lower action configuration with negative Boltzman 
factor which appears in the partition function is the spatial square of side 2 which has area Aj^a = 4. 
This shows that the positive and negative Boltzman weights are not balanced and that techniques 
such that the histogram method ||2^ can be applied. Equipped with the geometrical insight 
provided by this gauge invariant representation, we are working to design a suitable algorithm 
for simulating the loop fermionic action. The simplest case is QED in (1+1) dimensions or the 
Schwinger model for which all the Boltzman factors are positive i.e. one does not have to worry 
about computing the sign ah{S) . In addition, in a two space-time lattice there is a one to one 
correspondence between fermionic loops and the surfaces with border on them. So, to evaluate 
the path integral, the procedure is to generate surfaces with self-avoiding frontiers and with the 
constraint that they produce on each time slice open paths with ends of opposite parity. This 
was done recently by means of a Metropolis Monte Carlo algorithm and the results are very 
encouraging. 

Finally, let us mention that in this paper we only considered the simpler Abelian massless 
theory. In the non- Abelian case (see ref. for the loop formulation of the pure gauge theory), 
where there are different colors, the Pauli exclusion principle implies that the maximum number 
of such pairs at any site can not exceed the total number of degrees of freedom of the quarks. 
Therefore, more complicated diagrams arise. For instance, intersecting fermionic loops at vertices 
where there are more than a qq pair. Here the path to be followed is not unique. It is interesting 
to consider the massive case too. In this case, in addition to the dimers, we would have monomers 
produced by the mass term. The construction of a lattice path integral in terms of loops for full 
QCD is a task that will be considered in a future work. 
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APPENDIX I 

A general configuration J- with ajr = —1 obviously must include at least one connected 
fermionic loop contributing with a —1 to crjr. We will prove here that such connected fermionic 
path necessarily encloses an odd number of sites and thus, in D — 2, it is not possible to "fill" its 
interior region in such a way that all the sites are "visited" ( the condition required for a non-null 
contribution to Zp). 

Notice that, according to equation (p2[), the crjr associated to a configuration with a unique 
connected fermionic loop is negative if Ljrc /2 -I- Ajrc is an even number. 

In order to prove that all the fermionic loops of J-'^ which give a — 1 contribution to CTjf enclose 
an odd number of sites we need to express the number of its inner vertices as a function of its 
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length and area. In D = 2, a generic connected fermionic loop T'^ enclosing an arbitrary number of 
vertices can be obtained by a constructive procces begining with a plaquette p and then generating 
all the others fermionic loops by "appending" plaquettes to this diagram on each of its links in the 
two different ways showed in FIG. 6. 




(I) (11) 

FIG. 6. 



The variations of Njri 
as follows: 



and Ajr<z for the two cases illustrated in FIG. 6 can be summarized 



I) AA^c 

II) AA^c 



= 1 AL^e = 2, 

= 1 ALjpc = 0. 



From the above relations it is easy to check 



AA^c = 
AL^< 



(43) 



(44) 



where N^^ and A^p^^ denote the number of plaquettes of type / and // respectively. It is also easy 



to see that the number of inner vertices /y^ to a connected loop is equal to iVp^j- . 
Ajrc = 1 + AAjrc and Ljrc = 4 + Aijrc, 



2 



Therefore, since 



(45) 



From equations (g2|) and (^5|) we conclude that the contribution of Tf. to crjr, ajr^ is —1 if and 
only if ljr<: is an odd number. 

The same procedure extends to higher dimensions. (^) is true if the surface is built through 
steps of type /, // and /// where /// means adding a plaquette making contact in three links 
of the border. If the surface intersects itself, inner vertices belonging to the intersection lines are 
counted with the corresponding multiplicity. 



APPENDIX II 
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A. Structure of surfaces 



Consider a surface S contributing to the partition function. Its border is C = dS, a simple 
loop without intersections (not necessarily connected). This surface carries a sign: 

a{dS) = (-1)^-*+^ (46) 

where: 

N = Number of connected components in C 
L = Number of links in C 
A = Number of plaquettcs in S 

The sign of S depends only on OS since two surfaces with the same border differ by an even number 
of plaquettes. 

We divide S in spatial and temporal plaquettes: S = {St, Pt} where 

St = Spatial plaquettes of S in the section t 
Pt = Spatial links in the section t corresponding 
to temporal plaquettes of S between t and t+1 

Now wc have at each spatial section the paths dSt- Part of these paths propagates to t + 1, 
part comes from t — 1 and part does not propagate and belongs to C. 

Bt = Part of dSt that does not propagate 
= Part of dSt that propagates to t + 1 
= Part of dSt that comes from t — 1 



dSt = Bt ■ Rt ■ RT (47) 

Let us localize now the links of C. Its temporals links correspod to dPt- Its espatial links 
belong to cspatials plaquettes, in which case they belong to At, or to temporal plaquettes that can 
move up or down. 

— Spatial part of C in section t that propagates to i + 1 
= Spatial part of C in section t that comes from t — 1 

To keep track of all links in S (those that belong to any plaquette in S) we consider those 
joining two spatial plaquettes or two temporal plaquettes. The first ones are irrelevant to the sign. 
The second ones can join two temporal plaquettes in the same section (irrelevant) or in correlative 
sections: 

= Spatial links of S in section t that propagate^ to t + 1 

and come from t — 1 

To sum up (keys denote the set for all t and parenteses for a fixed t) 
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S^iSuPt} 
dSt = BtR+R^ 
Ct = {Bt,K+,K-) 

C^{Ct,dPt} 



Note that links in Bt, , are incompatible among them while lit, Rt ^ can share parts 
among them and with the former. 

There is a balance of creation and annihilation 

^(lA^I + |i?+|) = ^(|A7| + \R^\). (48) 
t t 

And what leaves at t — 1 is what arribes to t 

Pt-i - (nt_i, X+ 1, i?+ i) = {TluK^.R^). (49) 

It also verifies 



B. Localization of the sign 

The size of an objet X is |^| = number of links or plaquettes (depending if it is a path or a 
surface) taking into account its multiplicity. A path is fermionic if it consists of closed single lines 
and open single lines from even sites to odd sites. For these paths, P — {Pa^Qp} where a labels 
the open components and /3 the closed ones, we define the function 

(-i)^n(-i)^ (51) 

a 13 

The meaning of this sign is that it is the sign we get when we create a single line of lengh L using 
only the 'te operators. There are {L — l)/2 joining negative actions if the line is open and L/2 — 1 
if the line is closed. 

The sign in (28) is defined 

'^hiS) - ll{-l)^'iy{Pt-i ■ dStHPt) (52) 
t 

and it can be written 

ah{S)^Y[<j\Pt-,,PuSt) (53) 
t 

a\Pt-i,PuSt) = {-l)^^HPt-i-dSMPt) (54) 
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